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Abstract 

The future detection of gravitational wave forces us to consider the many ways in 
which astrophysics, gravitational wave theory and fundamental theory will interact. 
In this paper, I summarize some recent work done to develop such an interface. In 
' particular, I concentrate on how non-vacuum astrophysical environments can mod- 

, ify the gravitational wave signal emitted by compact binary inspirals, and whether 

■ signatures from the former are detectable by current and future gravitational wave 

^ ' detectors. I also describe the interface between gravitational wave modeling and 

, fundamental theory, focusing on the status of the parameterized post-Einsteinian 

framework (a general framework to detect deviations away from General Relativity 
in future gravitational wave data) and its current data analysis implementation. 



> 



1 Multi-Messenger Gravitational Wave Astrophysics 

I . The future detection of gravitational waves opens fertile ground to develop a collaboration between astro- 

^jq' physics, gravitational wave theory and data analysis. The traditional way in which such interdisciplinary 

work was thought to happen was fairly sequential: astrophysics would predict which sources gravitational 
wave theorist should model; theorist would provide highly accurate models for such systems; and data 
analysts would take such models and implement them in pipelines to maximize the probability of detect- 
' ing a signal. An accurate astrophysical system and waveform model are crucial as the first detections 

OS . will consist of signals that are completely buried in instrumental noise. 

^ ■ Once a signal has been detected, feedback in this cycle can occur. That is, data analysts can extract 

the best-fit parameters for the signal detected and, given enough detections, create population models 
O^l ' of the distribution of system parameters in the low-redshift universe. For example, given a sufficient 

number of compact binary signals, one could construct the distribution of masses, mass ratios, spins, etc. 
This information could then be returned to the astrophysics community to see how it matches predicted 
population models. 

Another route in which feedback can occur is through the data placing constraints on General Rel- 
ativity. Given a single detection, data analysts can determine the probability that the foundational 
\^ ' assumptions of General Relativity remain valid in strongly gravitating, highly dynamical systems, where 

tests of Einstein's theory currently do not exist. The most stringent tests we have today come from bi- 
nary pulsar observations (see e.g. [ij), but the gravitational compactness of such bodies, the characteristic 
mass of the system divided by its orbital separation, remains minuscule. Compact binary coalescences 
can in principle sample horizon-scale compactnesses, thus allowing for the most stringent tests yet. 

The interaction of gravitational wave theory with astrophysics is a two-way street. One way to proceed 
is to focus on how the dynamical evolution of a gravitational system affects the astrophysical environment. 
An example of this is the evolution of a circumbinary accretion disk when its binary black hole nucleus 
coalesces. In this process, the binary emits ~ 5 — 10% its rest mass in gravitational waves, depending on 
the orbital characteristics of the system. A change in rest mass leads to a modification in the Newtonian 
potential that the circumbinary accretion disk is bounded to, forcing the latter to readjust and leading 
to shocks in the material that then emits electromagnetic radiation (see e.g. [2I, [§]). 

A different way to proceed is to focus on how the astrophysical environment modifies the gravitational 
wave evolution. Let us revisit the example given above, ie. that of a circumbinary accretion disk with a 
comparable-mass binary black hole in a quasi-circular inspiral. One could imagine that the mere presence 



X 



^ Email address: nyuncs@physics.montana.edu 



2 



Gravitational Waves from Compact Binaries as Probes of the Universe 



of the accretion disk might modify the evolution of the binary, for example due to angular momentum 
transport by the disk or the disk's self-gravity. Surely, these and other effects will have an impact on the 
binary's inspiral rate, but are such effects measurable by future gravitational wave detectors? For the 
current example, the answer is likely no. That is, the effect of the accretion disk is probably too small to 
be detectable by LISA, because the comparable mass binary's inertia is too great relative to the effect of 
the disk. 

The effect of the astrophysical environment on the gravitational wave signal is not always negligible, 
however, and the prototypical system where the former matters is extreme mass-ratio inspirals. These 
systems consist of a solar-mass compact object spiraling into a supermassive black hole in a generic 
orbit. Due to the extreme mass-ratio, any small perturbation produced by a non-vacuum, astrophysical 
environment can have drastic effects on the small object's trajectory. Moreover, such systems can remain 
in the sensitivity band of low-frequency gravitational wave detectors for tens of years, depositing hundreds 
of thousands of radians of phase information in any given year. Thus, any small perturbation can 
accumulate over a year, leading to an observable effect. In this paper, I will describe a few papers 0-11] that 
have studied two possible astrophysical environment effects on extreme-mass ratio inspiral gravitational 
waves: the effect of accretion disks and the effect of massive perturbers. 

The interaction of gravitational wave theory and data analysis with fundamental theory can also pro- 
ceed in two different ways. Until very recently, the standard approach was top-down: one would pick a 
modified gravity theory, derive its equations of motion and solve them to find its predicted gravitational 
wave radiation. Once a prediction had been made, one could construct a template bank for this modified 
theory and then, given a gravitational wave detection consistent with General Relativity, constrain how 
large the parameters of the modified theory can be. Examples of this approach exist for Brans-Dicke 
theory 
gation 

violations of Local Position Invariance [2]| and the existence of extra-dimensions j22j • 

The problem with any of these approaches is that the construction of template banks and subsequent 
analysis of detected signals is computationally expensive for each alternative theory. Moreover, it is 
unclear which modified theory effect to concentrate on, since none of the above models is particularly 
more theoretically appealing than General Relativity. For example, in the massive graviton case, there 
is no self-consistent theory at the level of the action that would lead to such massive graviton effects. 

Lacking such a particular compelling model, a new framework to constrain deviations from Gen- 
eral Relativity was recently proposed: the parameterized post-Einsteinian (ppE) framework (23l |. This 
approach proposes to enhance the General Relativity template banks (spanned by system parameters) 
through the addition of new, well-motivated theory parameters. These parameters are such that when 
they acquire certain numerical values, they lead to waveforms predicted by all the modified gravity the- 
ories described above. The advantage of such an approach is that one need only concentrate on the 
construction of a single template bank. Matched filtering with elements from this bank would then allow 
the data to select the theory that best matches the signal, thus lifting the fundamental theoretical bias 
that General Relativity must be a priori correct. In this paper, I will summarize the current status of 
the ppE framework. 



l- fl2l. []3 , dyna mical Chern-Simons gravity 1^ lif. phenomenological massive graviton propa 
I, y, |l3|, [iTl |l8|, gravitational Lorentz- violation [l9|, and gravitational parity violation [20| 



2 Gravitational Wave Modeling 

Gravitational wave modeling has a long history and it would be quite pretentious and inappropriate to 
attempt a complete description here. Instead, I will only summarize the bits and pieces of gravitational 
wave modeling that will be necessary to understand the modeling of astrophysical and modified gravity 
effects, focusing on compact binary coalescences (inspiral-phase only) and ignoring any tidal deformation 
of the compact objects. I refer the interesting reader to the following review papers [24,] for post- 
Newtonian, comparable-mass inspiral modeling, |25j for self-force, extreme-mass ratio modeling, [26| 
for approximate, extreme-mass ratio modeling, and references therein. 
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2.1 Comparable-Mass Compact Binary Quasi-Circular Inspirals 



Let us begin with the modehng of gravitational waves emitted by comparable-mass, compact binary 
inspirals. The observed quantity at the detector is not the time-domain waveform, but the response func- 
tion: a quantity constructed from the contraction of the response tensor with the (transverse-traceless) 
gravitational wave metric perturbation: 



h{t) = F+{0s , 0s , 4^3)1^+ + Fx {es , <^s , ^s)/ix , 



(1) 



where /i+,x are the plus- and cross-polarized gravitational wave metric perturbations, while -F+_x are 
beam-pattern functions . These functions characterize the response of the detector to an incoming 
gravitational wave and are slowly-varying. 

The Fourier transform of the response function can be evaluated analytically if we assume that the 
phase varies on a much shorter time-scale than the amplitude, ie. in the stationary phase approxima- 
tion [23] . Let us define the cosine and sine responses for the dominant £ = 2 mode 

h{t) = hc{t) + hs{t) , hc{t)^ A Qc{i,P) cos2(i3, hs{t)^AQsii,l3)sin2(l,, (2) 

where A = —M./Dl {2nA4F)'^^'^ is a slowly- varying amplitude, while 

Qc{l, /3)=2{1 + cj) C2pF+ - 4c,S2pF^ , Qsih /?) 2 (l + c?) s2pF+ + 4c,C2^Fx , (3) 

and Ci — cos t, C2p = cos 2/3 and = sin 2/?, with t and /? the inclination and polarization angles 
respectively. Then, the Fourier transform in the stationary phase approximation is 



hcif) = - 



1/2 



'2/3 



Qc{i,P) [2F{to) 



-7/6 ^-i(*+7r/4) 



^384/ Dl 

while hsif) — ihc^sif)- The phase of the gravitational wave Fourier transform is 

^F{to)]=2n J 



(4) 



(5) 



The quantity to is the stationary-point, defined by 20(to) = 27r/, essentially the region in which the 
Fourier transform does not average out. In Eq. is the luminosity distance, M = rf^^m is the 

chirp mass, with 77 — mim2 / m'^ the symmetric mass ratio and m = mi + m2 the total mass, while F is 
the orbital frequency. 

Gravitational wave detectors are the most sensitive to the phase of the response function, Eq. ([5]), 
whose functional form is controlled by the orbital frequency's evolution equation. In vacuum General 
Relativity, the orbital frequency can only evolve due to the emission of gravitational waves. Using the 
balance law E — — {Cqw), i-e. the rate of change of the binary's binding energy is exactly balanced 
by the flux of gravitational wave energy-momentum suitably averaged over several wavelength, one can 
derive 

to leading order in 2ttA4F ~ w'^, where v is the binary's orbital velocity. Of course, if Kepler's third 
law is modified or if additional sinks of energy- momentum are present, for example due to scalar-field 
emission or to accretion-disk induced angular momentum transport, then the gravitational wave phase 
will be modified accordingly. Ignoring these effects for the moment, the Fourier transform is then 



^(circ) - (^334) 



1/2 



-2/3 



X5/6 

Dl 



-7/6 p 



exp 



2nft, 



4 



128 



inMf) 



-5/3 



(7) 



where Q — Qc + iQs and we have retained only the Newtonian term in the phase. 
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2.2 Extreme Mass-ratio, Compact Binary Quasi-Circular Inspirals 

Let us now consider extreme-mass ratio inspirals in a quasi-circular orbit confined to the equatorial plane 
of the supermassive black hole. The dominant ^ = 2, gravitational wave mode accumulates a phase 



= 2 



m dt = 2 [\r)^ = lMsMBH^5/2 
A/-T„,, Jro ^ 16 mco ^ 



(8) 



where Tobs is the observation time, 51 is the orbital angular frequency, -Msmbh is the supermassive black 
hole mass, mco is the stellar-mass compact object's mass and {Rq, Rf) are the AfsMBH-normalizcd initial 
and final orbital radius respectively. The quantity t here is the dimensionless observation time defined 
by T = (256/5) (mco7obs/.^^sMBH)- ^^"^ '"^^ verify that for low-frequency gravitational wave detectors, 
with peak sensitivity in the micro- to deci-Hz band, the accumulated gravitational wave phase can reach 
millions of radians for separations in (i?o,i?/) € (30Msmbh, ^isco)j where i?isco is the radius of the 
innermost stable circular orbit for a test-particle on a geodesic of the background spacetime. 

The accumulation of gravitational wave phase information clearly depends on the rate of inspiral, f. 
This quantity can be computed from r = E/{dE/dr), which again depends both on the binary's binding 
energy E{r) and its rate of change. If additional sinks of energy are present in the problem or if the 
binding energy is deformed, then the accumulated gravitational wave phase will be different than what 
one would expect for a vacuum inspiral. 

The accumulated gravitational wave phase in Eq. ([5]) can be used as a quick measure of how im- 
portant non-vacuum or non-General Relativity effects are, but this measure is only approximate since 
it uses Newtonian relations, which are highly inaccurate for extreme-mass ratio inspiral. A much more 
appropriate scheme is that of the effective- one-body approach, recently implemented for extreme-mass 
ratio inspirals in [28l, ^295 and references therein. In this scheme, one maps the two-body inspiral into an 
effective inspiral of a test-particle with mass rj in orbit around a supermassive black hole with mass to. 
The equations of motion then reduce to the Hamilton-Jacobi equations 

dr _ A{r) dH'"''^ d^ _ dH""^^ 

'dt ^ ^/D{r^ dpr, ' 'dt ~ dp^ ' 

dpr, _ A{r) ai/™^^ dp^ 



dt -s/lW) ' 



= 7"$ . (10) 



where jj™^' = H'^°'^^/fi is the reduced (i.e., dimensionless) real Hamiltonian and Tip/^ is a reduced 

radiation-reaction force that controls the rate of inspiral (there is no radial component on average here 
because we are dealing with an equatorial, quasi-circular orbit). The azimuthal radiation- reaction force 
can be modeled via J^j — rj^^L = —rf^^LQ^f^ — — (r/J7)^^i?GW5 where L is the rate of change of the 
binary system's orbital angular momentum, while Xqw and -Egw a-re the angular momentum and energy 
radiated in gravitational waves. 

The energy flux is composed of several contributions. The most accurate model for the average, energy 
fiux radiated to infinity is the factorized resummation 

i=a ■m=e 

EGW^-ry2y2{rnn)'\Rhi^\\ (11) 
Idtt ■^-^ 

1=1 m=l 

where R is the distance to the observer and hem is the harmonically-decomposed gravitational wave metric 
perturbation. In addition to this, there is also energy flux into any trapped surfaces, black hole horizons, 
that can be modeled by solving the Teukolsky equation perturbatively in a post-Newtonian expansion 
(see e.g. (s^llll). Of course, if there are additional sinks of energy, these must also be accounted for here 



by adding contributions to the total energy flux. 
2.3 Gravitational Wave Signatures 

In the previous sections, we have briefly summarized some salient features of gravitational wave mod- 
eling for comparable-mass and extreme mass-ratio quasi-circular inspirals, but are modifications in this 
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modeling observable? A plausibly criterion to decide whether a certain gravitational wave modification 
is detectable by a given detector is the following: 



> 



10/ p ifp>10, 







if p < 10 , 



where is the square of the signal-to- noise ratio (SNR), defined as 



p\h)^4 



df 



Shif)' 



(12) 



(13) 



As before, h is the Fourier transform of the measured gravitational wave response function, while Sh{f) 
is the spectral noise density curve of the detector. The piece- wise nature of this criterion accounts for 
the fact that there is a threshold SNR below which the signal cannot be detected in the first place. 

The above requirements give us a rough sense of when a certain modification might be detectable, 
but can such a modification be distinguishable from a vacuum waveform. To address this question, one 
can consider the SNR of the waveform difference 



p (Sh) — min 

^2 



df 



SnW 



(14) 



where hi and /12 are the Fourier transforms of two waveforms (the "signal" and "template"), normalized 
such that p{hi) — p(/i2) = 1- The template depends on parameters A2 that may be different from the 
true astrophysical ones, where the minimum difference corresponds to the best fit. One can then consider 
this SNR, minimized over template parameters, to construct a measure similar to that of Eq. (jl2p . 
The minimization will partially account for the effect of possible template parameter degeneracies with 
non- vacuum or non-General Relativity effects. 



3 Astrophysical Imprints 

In this section, we consider two astrophysical environment effects in the modeling of gravitational waves. 
We begin with massive perturbers, summarizing the results of [4], and then move on to accretion disk 
effects, summarizing the results of |5|, |6|. I refer the interested reader to these papers and the many 
references therein, which will not be included in this paper for simplicity. 

3.1 Massive Perturbers 

Let us consider an extreme-mass ratio inspiral, where the stellar-mass compact object is in a quasi- 
circular, equatorial orbit around a spinning supermassive black hole. Let us now imagine that there is 
a third compact object, with mass Mgcc and at a distance rsec from the supermassive black hole. We 
here imagine that Msoc is comparable to Msmbh and much greater than mco • This scenario is depicted 
in Fig. [TJ Clearly, this is a three-body system, composed of two sub-binaries: the extreme-mass ratio 
system and the Msmbh — -Mgcc system. Because of this, the extreme mass-ratio gravitational waves will 
be emitted from an accelerated frame (a frame that rotates due to the presence of the secondary). 

The modifications to the trajectories of the small body can be modeled in the framework of the 
effective-one-body scheme described in Sec. 12.21 The dominant effect is simply a Doppler shift in the 
extreme mass-ratio inspiral gravitational wave frequencies, which then leads to an integrated modification 
in the gravitational wave phase. The implementation of this correction is simple: divide the right-hand- 
side of Eq. (|TU| by the appropriate Doppler factor 

$ = f7 ^ $ = [l + uios(t,5==7r/2)] , (15) 

where wios is the velocity along the line of sight 

wios(i) = ( ) WNowtCOs(cjNowti + <^) sin(i) . (16) 

V -'WTot / 
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Figure 1: Schematic view of the extreme- mass ratio system (in the xy plane), the massive perturber Mscc 
(at a distance rgcc), and the line of sight, t is the inclination between the primary-secondary supermassive 
black hole's orbital angular momentum vector and the line of sight. 



with WNowt — {GM^ot/rscc)^^^ the Newtonian virial velocity, r^Ncwt = (CMTot/^'scc)^^^ the Newtonian 
angular velocity for an object in a circular orbit, and 6 an initial phase offset, with (f2Ncwti + ^) the orbital 
phase of the Msmbh — Afsoc system. A constant relative speed can be re-absorbed in a redefinition of 
the masses, and thus, it does not lead to a measurable correction. In Eq. we have ignored the 

appropriate Lorentz factor F, since i^Newt/c 1, and we have removed the constant velocity drift in vios 
by choosing S = Tr/2. 

This simple modeling neglects other gravitational wave generation effects. One could incorporate such 
effects by introducing an external, vectorial force (the product of the the total mass of the system and 
the time derivative of the velocity of Eq. (IT61) ') to Hamilton's equations. Such a modeling would require 
a non-adiabatic evolution, which would lead to eccentric and inclined extreme-mass ratio trajectories. 
The magnitude of this correction, however, is small, as it scales with the tidal force of the perturber 
on the center of mass of the extreme mass-ratio inspiral, relative to the stellar-mass compact object's 
acceleration due to the secondary supermassive black hole. Since the tidal force scales as ~ Msec/''seci 
this effect is suppressed relative to the acceleration by a factor of r/rsec ~ 10~^ for an extreme- mass ratio 
inspiral with orbital separation SOMsmbh and a primary-secondary binary orbital separation of 0.01 pc. 
This ratio would not be small if the perturber were in the near-zone (less than a gravitational wavelength 
away) from the center of mass of the extreme-mass ratio system. 

We can now estimate what the distance and mass of the secondary black hole must be for its ef- 
fect on extreme mass-ratio inspiral gravitational to be detectable by a LISA-like mission. The time- 
dependent Doppler shift described above leads to a correction in the gravitational wave phase of A0gw = 
ViogToitsN / {2c) ^ where N is the number of radians in the waveform, Tobs is the observation time. Let us 
further define e as the detectable fractional phase shift: e = A0gw, detect /A^, a fiducial value of which is 
e = 10~^, or 0.1 radians over ^ 10^ radians for a typical one- year inspiral. 

The definition of e is a curve in rgec ~ -^^Sec space that delineates the secondary parameters that 
would lead to a sufficiently large dephasing to be detectable. This is shown in Fig. [2] for a range of e 
values. Any secondary system with mass and separation above the curves would be observable. The solid 
curves plot the relation described in the previous paragraph, while the dashed curves account for the next 
v\ost'^ order effect. For comparison, the region of (Mscc, ^Sec) space that falls in the pulsar-timing-array 
sensitivity band are also shown. Of course, the gravitational waves that would be detectable by pulsar- 
timing-arrays would be those generated by the Msmbh — Mgcc binary, and their distance to Earth would 
have to be sufficiently small. In principle, however, this scenario allows for the possibility of coincident 
future detection of gravitational waves from LISA-like detectors and pulsar-timing-arrays. 

Just because an effect is sufficiently large to be in principle measurable does not imply that it can 
be distinguished from other vacuum effects. One must therefore worry about possible degeneracies of 
vacuum system parameters with those introduced by the secondary perturber (Afsec and rgec in this 
case) . We can get a sense of whether this is the case by computing the Fourier transform of the modified 
gravitational wave in the stationary-phase approximation. Following Sec. 12.11 one finds the following 
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Figure 2: Range of secondary masses and separations that could be measurable by LISA given a suffi- 
ciently strong EMRI. The region above the solid and dashed lines would be observable. Measurement 
of the leading-order effect gives a determination of the combination (MsccSint)/rg(,^, while measuring 
the next-order effect gives a determination of the combination (Mgoc sin t)'^/^/r'g^^. Thus measuring both 
effects together allows both Msec sin l and rsec to be determined. 



correction to the phase 



A* 



-5/3 
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' Sec 



1094713344 



' Sec 



(17) 



where we have neglected a constant term that is no measurable. Observe that the corrections introduced 
by Msoc cannot be re-absorbed by a redefinition of vacuum system parameters, which implies that 
the effect of the secondary are (at worst) weakly-correlated vacuum terms. Moreover, if one could 
independently constrain the coefficients in from of the dominant f~^^^ term and also the /~^^/^ and 
terms, one would be able to break the degeneracy between Msec and rsec in the leading-order term. 



3.2 Accretion Disks 

Let us now consider a different astrophysical environment: an extreme mass-ratio system in a circular 
equatorial orbit that is embedded in the accretion disk associated with the supermassive black hole. We 
refer the interested reader to for astrophysical scenarios that could lead to this configuration. Of 
course, the effect of such a disk on the extreme-mass ratio inspiral will depend on the disk properties, 
which are not accurately known. One must then think of all possible effects that could be present and 
determine which one leads to the largest correction to extreme-mass ratio inspirals. 

The dominant effect is that induced by the gravitational torque induced by spiral density waves in 
the accretion disk (ie. migration torque in planetary dynamics). Migration can occur whether a gap is 
cleared out or not, with the former leading to the largest effects. One might naively think that the stellar- 
object's mass accretion would also greatly modify the object's motion, but this is not the case because of 
quenching by a variety of processes, predominantly limited gas supply and radiation pressure. Migration 
then becomes dominant, but its effect on gravitational waves still depends on the disk's properties. If one 
models the latter as an a (Shakura-Sunyaev) disk (where viscosity is proportional to the total pressure), 
then all disk effects become negligible. But if the disk is better described by a /3 profile (where viscosity 
is proportional to gas pressure only) in the radiation-dominated (extreme-mass ratio) regime, then the 
effect of migration would be measurable. Therefore, if such an accretion disk effect is detected one would 
know that the disk associated with that particular extreme mass-ratio binary is of /3 type. 

Figure |3] summarizes the impact of the main accretion disk effects on the gravitational wave. This 
figure plots the dephasing (the difference between the gravitational wave phase computed with and 
without disk effects) as a function of the final orbital radius, assuming a one year inspiral. The left 
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and right panels show results for different extreme mass-ratios, while the different curve colors and 
styles correspond to different disk effects: black is for accretion onto the small compact object, blue is 
for migration and green is for azimuthal winds. The latter corresponds to a difference in the orbital 
compact object velocity and the gas velocity, which can push the former forward or backward. The 
total accumulated gravitational wave phase is denoted by a thin magenta line, while the thick magenta 
lines correspond to two (an optimistic and a pessimistic) measures of the accuracy of a classic LISA-like 
mission. Notice that migration can lead to dephasing of order 10'^ — 10"^ radians for /3 disks, while the 
same effect is essentially negligible for a disks. Therefore, a measurement of such disk effects could allow 
us ti 
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Figure 3: Gravitational wave dephasing as a function of final radius in units of Msmbh induced by 
different disk effects. Solid (dotted) curves correspond to a {(3) disks, with different colors indicating 
different disk effects: black corresponds to Bondi-Hoyle-Lyttleton (BHL) accretion onto the small object, 
green to azimuthal winds and blue to migration. The thin, solid magenta line is the total accumulated 
gravitational wave phase in vacuum. The thick, solid (dashed) magenta line corresponds to a measure 
of the accuracy to which LISA can measure the gravitational wave phase for a source at 1 Gpc (10 
Mpc). Observe that certain disk effects, like migration, can leave huge imprints on the gravitational wave 
observable, inside the LISA accuracy bucket, provided the disk is a /3 one. 

As before, a simple dephasing argument does not guarantee that disk effects can be disentangled 
from other vacuum system parameters. One can get a sense of whether such disentanglement is possible 
by computing the Fourier transform of the gravitational wave response function in the stationary-phase 
approximation: 

*/*vac = 1 - il«r™SMBH,l^^l^BH,5C<' ' (18) 

and 

\h\/\hU - 1 - Bia'i^^^SMBH,lMll,BH,,d'4' : (19) 

where ai = a/0.1 is the normalized a-viscosity parameter, tosmbh,i = "t-smbh/0.1 is the normalized 
supermassive black hole accretion rate rh = AfsMBH/-^SMBH,Edd with Msmbh, Edd the Eddington rate. 
Similarly, AfsMBH,5 = A^smbh/(1O^M0) is a normalized supermassive black hole mass, qq = is 
the normalized mass-ratio q = Msmbh/toco and uq = (7rA^/)/(6.15 x 10~^) is a normalized reduced 
frequency for a gravitational wave frequency of 10~^ Hz. The parameters (Ai, Bi,di,Ci) are given in 
Tabled] Of course, these expressions are only valid when the accretion disk effects are small perturbations 
away from the vacuum evolution (ie. at sufficiently small separations). 

The important result here is that the dependence of the disk effects with frequency, ie. the parameter 
as, is of opposite sign with respect to those that arise in the vacuum gravitational wave phase. This is 
because disk corrections are largest for large radii, equivalent to weak-field vacuum effects, while post- 
Newtonian corrections to the vacuum gravitational wave phase are largest for small radii, equivalent 
to strong-field effects. This suggests that migration effects are weakly correlated with vacuum system 
parameters, as one could not reabsorb these disk effects by choosing different vacuum parameters. 
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^1 


Bi 


as 




05 


Cl 


C2 


BHL a 


3 (-8) 


2 (-7) 


1 


4 


-20/3 


-1 


-5 


BHL (3 


l(-5) 


l(-4) 


6/5 


79/25 


-79/15 


-4/5 


-17/5 


W a 


6 (-17) 


1(-16) 


1 


16/5 


-16/3 


-1 


-3 


w p 


6 (-12) 


4(-ll) 


6/5 


59/25 


-59/15 


-4/5 


-7/5 


Ml a 


3 (-10) 


4 (-9) 


1 


16/5 


-16/3 


-1 


-3 


Ml P 


l(-6) 


3 (-6) 


6/5 


59/25 


-59/15 


-4/5 


-7/5 


M2 asc 


8 (-6) 


2 (-5) 


1 


-2/5 


-8/3 





1 


M2 ^sc 


6 (-3) 


2 (-2) 


1/4 


-1/8 


-25/12 


1/2 


5/8 


M2 /3ipp 


6 (-4) 


2 (-3) 


4/7 


-17/70 


-7/3 


2/7 


11/4 



Table 1: Columns are parameters in Eqs. (jl8p and (fT9)) and rows are disk effects. The notation x [y) = 
a; X lO*' in radians for Ai and dimensionless for Bi. Observe that the frequency exponent 05 < 0, implying 
that these accretion disk effects are dominant at small frequencies (large radii). 



4 Modified Gravity Imprints 

In this section, we discuss how modified gravity imprints can manifest in the gravitational wave observable. 
We describe a scheme developed to constrain such imprints, the so-called parameterized post-Einsteinian 
(ppE) framework. We mainly follow the presentation in [16| and the recent results of [32|. 



4.1 Top-Down Approach 

One approach to testing Einstein's theory with gravitational wave observations is a top-down one: pick a 
modified gravity theory, study its consequences and match-filter gravitational wave data with a template 
bank specifically constructed for that theory to test for possible deviations consistent with that theory. 
Most studies carried out so-far considered only quasi-circular coalescences in the inspiral phase. For such 
waveforms, the modified Fourier transform of the response function in the stationary-phase approximation 
can be modeled as 

M/) =/^GR(/)(l + «^^")e'''"\ (20) 

where the GR transform was already described in Sec. 12.11 while (a, 6) are numbers and (a,/?) are 
functions of system parameters and parameters of the theory. We present their values in the Table [2 



Theory 
Brans-Dicke 
Massive Graviton 
Parity Violation 
Non-Dynamical CS Gravity 
G{t) Theory 
Lorentz Violation 

Conservative EDGB 
Extra Dimensions 






^v5A 
AnS-dcs^A 





EDGB 







1 

1 

_8 
3 



/3 

._5 SL„2/5~ 

3584 WBD " 






75 



„ 65536' 

25^-4/5a 

raV Cedgb 



,__iM^-4(3_26,?_ 



2477^) 



13 

■ 3 



-aLv - 1 



3 

13 



Table 2: Parameters that define the deformation of the response function in a variety of modified gravity 
theories. The notation • means that a value for this parameter is irrelevant, as its amplitude is zero. All 
other parameters are defined in the text. 



Let us explain how the symbols of Table [5] are defined. In Brans-Dicke theory, 5* is the difference in the 
square of the sensitivities and cjbd is the Brans-Dicke coupling parameter In the phenomenological 
massive graviton theory, D is a certain distance measure, z is redshift and \g is the Compton wavelength 
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of the graviton In the parity violation case, the amplitude deformation depends on the inclination 
angle and the beam pattern angles via (l^ 

(F^+F^)c.(l + cf) 



while D is a measure of parity violation. In non-dynamical Chern-Simons (CS) gravity [33| , the quantity 
S-dcs is the time derivative of the CS scalar field (ie. the magnitude of the canonical embedding coor- 
dinate) In a varyiiig Newton's constant theory, Gc is the value of the time derivative of Newton's 
constant at coalescence [2l| . In the phenomenological Lorentz- violating theory, Alv is a distance scale at 
which Lorentz- violation becomes important, while 7 is the graviton momentum exponent in the defor- 
mation of the dispersion relation [l9|. In Einstein-Dilaton-Gauss-Bonnct gravity, Cedgb is the coupling 
parameter in the theory, and we here present only the modification in the response function introduced by 
the deformation of the background Schwarzschild metric [s^. In theories with extra dimensions, dM/dt 
is the mass loss due to gravitons leaking into the extra dimension. 

As one can see from the above table, there are many possible deviations from General Relativity 
that one could consider, and none of the above has been sufficiently well-studied to determine whether 
any of them are compelling alternatives. A top-down approach would require us to construct as many 
template banks as there are alternative theories to consider. Of course, each of these banks would contain 
additional dimensions in parameter space (and thus many more templates), as one would have to allow for 
variations of the fundamental parameters of the modified theories. This might increase the computational 
cost beyond current capabilities. 

But perhaps the biggest drawback of a top-down approach is that one must pick a theory before- 
hand, and this does not allow the data to select the theory that fits it best. Instead, this choice is made 
by the gravitational wave theorist or data analyst. It is entirely possible that, if General Relativity is 
incorrect in the strong field, none of the above theories will represent the correct modification. Thus, if 
one implements a top-down approach, one is forcing a certain fundamental bias into the analysis, ie. the 
fundamental bias that we, as theorists, know what modifications of gravity are possible. 

4.2 Bottom-Up Approach 

A new scheme has been proposed to alleviate this fundamental bias: the ppE framework [l6| . In this 
scheme, one enhances General Relativity templates /igRi which depend on system parameters Agys, 
through the introduction of theory parameters AppE: 

h{f; Xsys) ^ hGKif; Xsys) + Shif; Asys; AppE) . (22) 

When these parameters acquire certain values, the templates reduce exactly to those predicted in Gen- 
eral Relativity, ie. Sh{f;Xsys',^p^) = 0. When they acquire other values, the templates reduce to the 
predictions of modified gravity theories. The idea then is to match filter or perform Bayesian statistics 
with this enhanced template bank to allow the data to select the best-fit values of AppE that best fit the 
signal. 

This framework emerges from a generalization of the parameterized post-Newtonian scheme, developed 
in the 1970's to test modified gravity theories with Solar System experiments 35-4^. That scheme 



proposed the enhancement of the main quantity which all Solar System observables depend on, the 
weak-field expansion of the metric tensor in a certain coordinate system and gauge, by theory (ppN) 
parameters: 

9tiu{x^] Asys) gtiu{x^; Asys, AppN) • (23) 

When these parameters vanish, then the metric tensor reduces exactly to that predicted by General 
Relativity g^,y(a;^; Asys, 0), while when they don't, they describe the weak-field expansion of the metric 
tensor in a plethora of modified gravity theories. The idea is then to constrain the values of AppN by 
measuring Solar-System observables that depend on the metric tensor. 

The ppN and ppE frameworks, however, differ in the priorities that drive the parameterized General 
Relativity deviations. While the former does not care about how many theory ppN parameters are 
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introduced, in the latter we wish to minimize the number of ppE parameters. Otherwise, the inclusion 
of too many ppE parameters dilutes the information in the extracted best-fit, raising the false-alarm 
probability. One must then find a balance between the number of parameters introduced and the amount 
of bias contained in the templates. 

The precise form of ppE templates clearly depends on the system under consideration and the level 
of theoretical complexity that one wishes to include. Keeping in line with the rest of this paper, we 
concentrate here on non-spinning, quasi-circular inspirals. In their original work, 16!] also proposed ppE 
templates that contain the merger and ringdown phases in a hybrid fashion, but we will not discuss these 
here. As for the theoretical complexity, this depends on the number of additional ppE theory parameters 
one wishes to include. 



The simplest ppE realization, which we refer to as the restricted ppE templates , is the following [41ll4 



,5M/;Asys;a,/3) = /iGR(a^^°™)e"'"''™ , (24) 

where the theory parameters (a,/?) are to be searched over, while (apN,^PN) are numbers equal to 
some of the post-Newtonian predictions for amplitude and phase frequency exponents, ie. (apN,&PN) 
are not allowed to vary. For example, a Newtonian prediction for the amplitude and phase frequency 
exponents would be aNcwt = —7/3 and &Nowt = —5/3, while a 1.5PN prediction would be ai,5PN — 1 and 
bi.5PN = —2/3. 

One can think of these tests as consistency checks of the post-Newtonian expansion. That is, given 
a detection of a gravitational- wave, one can in principle measure the chirp mass and the mass ratio 
by detecting the Newtonian and IPN corrections to the phase. Allowing for a 1.5PN parameterized 
deviation in the phase then allows one to check whether this term is indeed as predicted by the post- 
Newtonian expansion. Alternatively, one can also think of these tests as straw-men indicators for GR 
consistency/inconsistency. That is, given a low signal-to-noise ratio detection, one can ask whether the 
signal detected is consistent with a General Relativity hypothesis (a, /3) = (0, 0) jisf . 

The great disadvantage of the restricted ppE formalism is that, although one might be able to discern 
whether deviations are present, one cannot tell what kind of deviations those are. To alleviate this 
problem, one can introduce free ppE templates , the simplest version of which is [l6j 

8h{f- Asys; a, /3) = hoK {a u'^) e^^"' , (25) 

where now (a, /3, a, b) are free ppE parameters. Observe that different choices of ppE parameters can 
recover all the modified theory predictions listed in Table [2] Furthermore, one can show that a generic 
power-series deformation of the binding energy or the gravitational wave luminosity will lead to a metric 



with the above ppE parameterization 16|. 



Given a stretch of gravitational wave data, one can then search via matched filtering with free ppE 
templates and allow the data to select the (a, /3, a, b) ppE parameters that fit it best. If a deviation is 
present, such a procedure would not only signal the presence of a deviation, but it would also signal what 
type of phase or amplitude deformation is introduced. That is, one would know the best-fit numerical 
values of a and 6, which encode what type of modified gravity effect is present. For example, as suggested 
by TablelU a non-zero detection of /3 associated with b = —7/3 would signal the presence of a scalar dipolar 
mode, while one associated with b = —13/3 would signal the existence of an anomalous acceleration. 

We have recently carried out a detailed Bayesian analysis with free ppE waveforms, whose main results 
are shown in Fig. |4]^32]. This figure shows the (3CT)-bounds on /3 for different values of b (ie. an exclusion 
plot for anything above the dotted curves), given gravitational wave detections consistent with General 
Relativity. The different dotted curves correspond to different binary systems, where we consider the 
inspiral phase only. For comparison, we also plot current bounds on (3 from the double binary pulsar, as 
well as a few Solar System constraints. Observe that for b > —2, gravitational wave observations could 
rule out a large sector of (/3, b) parameter space that is currently unconstrained. These results where 
confirmed using both a Fisher analysis as well as mapping out the likelihood surface with a Markov- 
Chain Monte-Carlo techniques. 

Given a gravitational wave detection, one can also ask whether a General Relativity waveform (an 
Einstein hypothesis) or a non-General Relativity waveform (a non-Einstein hypothesis) best fits the data. 
This can be established by computing the Bayes factor or the odds-ratio of these two hypothesis. When 
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Figure 4: Left: (3(T)-bounds on /3 for different values of b for a single SNR = 20 aLIGO/aVirgo detection. 
The two dotted lines correspond to two sources with different mass ratios, total masses, and sky locations, 
but both at redshift z — O.l {Dl — A62Mpc). The solid line is the (3(T)-bound on (3 from the golden 
pulsar (PSR J0737-3039), while other symbols are bounds from Solar System experiments and other 
aLIGO analyses. Right: Same as left but for two classic-LISA sources at redshift z = 1 and z = 3. 

using ppE templates, the hypothesis become nested, but the above described analysis can still be easily 
carried out. The Bayes factor will increase the larger (a, /?) are, and so assuming a threshold Bayes factor 
of 10 (1 to 10 betting odds that one of the two hypothesis is incorrect and the other is correct) one can 
derive the accuracy to which (a, /3) can be measured as a function of b. Implementing such a scheme, 
one obtains similar results to those already plotted in Fig. |4l 

5 A Roadmap for the Future 

The era of multi-messenger gravitational wave astrophysics is at our doorstep. A full exploitation of 
gravitational wave information will require a deep collaboration between astrophysicists, gravitational 
wave modelers, general relativists and data analysts. In this paper, I have described recent efforts to 
extract information about the astrophysical and theoretical environment in which compact binaries might 
be evolving. 

On the astrophysical side, we have seen that the presence of a secondary perturber or an accretion 
disk can lead to observable effects in extreme-mass ratio inspiral gravitational waves. As for the former, 
these effects become non- negligible when secondary perturbers are at a distance of 0.1 pc or less from 
the center of mass of the extreme- mass ratio binary, given secondary masses of IO^Mq. As for the 
latter, angular momentum transport due to density waves, analogous to planetary migration, can leave 
observable signatures in extreme-mass ratio binaries embedded in the thin disks, provided the latter are 
of /3-type. These effects seem weakly degenerate with other system parameters, as suggested by a simple 
Fourier analysis. 

On the fundamental theory side, we have seen that a top-down approach leads to a variety of wave- 
form predictions. There is no particular modified gravity theory that is more compelling than General 
Relativity, which forces us to consider all theories on an equal footing. Performing a data analysis study 
for all such theories might be computationally prohibitive. Instead, one can carry out a bottom-up ap- 
proach by match- filtering with ppE templates. These waveforms enhance General Relativity templates 
with certain well-motivated theory parameters, which then allows the data to select the parameters that 
best fit it, lifting some degree of fundamental bias. 

A variety of paths present themselves for future studies. One could, for example, consider other 
astrophysical environment effects, or other specific top-down modified gravity theories, to discover what 
the corrections to the waveform are. Alternatively, one could consider whether the presence of an as- 
trophysical environment could prohibit tests of Einstein's theory. This seems quite unlikely, given that 
modified gravity effects if present should be there always, while astrophysical environment effects should 
only affect a hand-full of sources. Finally, a realistic implementation of the ppE idea will require the gen- 
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eralization of the free ppE model to more complex and accurate waveforms. If this is not accomplished, 
mismodeling bias is likely to contaminate ppE studies, hindering the ability of gravitational waves to 
constrain Einstein's theory. Work along this lines is promising. 
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